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We use the coupled cluster method to high orders of approximation in order to calculate the 
ground-state phase diagram of the XXZ spin-s kagome antiferromagnet with easy-plane anisotropy, 
i.e. the anisotropy parameter A varies between A = 1 (isotropic Heisenberg model) and A = 0 
(XY model). We find that for the extreme quantum case s = 1/2 the ground state is magnetically 
disordered in the entire region 0 < A < 1. For s = 1 the ground state is disordered for 0.818 < A < 1, 
it exhibits %/3 x %/3 magnetic long-range order for 0.281 < A < 0.818, and q = 0 magnetic long- 
range order for 0 < A < 0.281. We confirm the recent result of Chernyshev and Zhitomirsky ( Phys. 
Rev. Lett. 113, 237202 (2014)) that the selection of the ground state by quantum fluctuations is 
different for small A (XY limit) and for A close to one (Heisenberg limit), i.e., q = 0 magnetic 
order is favored over \/2> x \/3 for 0 < A < A c and vice versa for A c < A < 1. We calculate A c as 
a function of the spin quantum number s. 

PACS numbers: 75.10.Jm, 75.10.Kt, 75.50.Ee, 75.45.+j 


Introduction. The investigation of the ground state 
(GS) of the quantum antiferromagnet on the kagome lat¬ 
tice is one of the most challenging problems in the field 
of frustrated quantum magnetism. Over many years nu¬ 
merous theoretical methods has been applied to under¬ 
stand the GS properties of the kagome antiferromagnet 
(KAFM), see, e.g., Refe.QHall and references therein. 

While it became clear very early that GS magnetic 
long-range order (LRO) is absent for the s = 1/2 Heisen¬ 
berg KAFM, there was a longstanding debate on the na¬ 
ture of the non-magnetic quantum GS. Recent large-scale 
numeric s 14 ’ 15 ' 19 provide strong arguments for a gapped 
Z2 topological spin-liquid GS for spin quantum number 
s = 1/2. 

Other recent investigations have been focused on 
higher spin s > 1/ 2 18 ' 24 ' 27 ~ — and also on the anisotropic 
XXZ-model . 10 ' 24 ' 25 Both modifications have relevance for 
the experimental research, see, e.g., Refs. HT^ - lTol More¬ 
over, anisotropic spin models are of great interest with 
respect to engineering models of quantum magnetism on 
optical lattices, see, e.g., Refs. l4ll-|43l Since higher spin 
quantum numbers as well as spin anisotropy, in general, 
lead to a reduction of quantum fluctuations, see, e.g., 
Refs. [~T~il fbil GS magnetic LRO for the KAFM might be 
facilitated. However, for the isotropic s = 1 Heisenberg 
KAFM there is clear evidence that there is no magnetic 
LROj 18 i 28 ~ — (Note, however, that in Ref. H3 v/3 x \/?> 
GS LRO for integer spin quantum numbers including 
s = 1 was reported.) For s > 1 several approaches 
lead to indications for >/3 x \/Z GS LRO.^ ' 18 ' 24 ' 27 On 
the other hand, recent density matrix group (DMRG) 
calculations^ have demonstrated that for s = 1/2 the 
XXZ KAFM remains in a magnetically disordered GS 
for the entire range of the anisotropy parameter A be¬ 
tween the XY point (A = 0) and the isotropic Heisenberg 
point (A = 1). 

Motivated by the recent Letter of Chernyshev and 
Zhitomirsky^! we use the coupled cluster method (CCM) 


in high orders of approximation to calculate the s — A GS 
phase diagram of the spin-s XXZ KAFM with easy-plane 
anisotropy. The corresponding Hamiltonian is 

H = X ( s i s j + s i sV j + As i s l) . 0 < A < 1, (1) 

[i,j] 

where the sum runs over all nearest-neighbor pairs. The 
CCM is a very general ab initio many-body technique 
that has been successfully applied to stro ngly frustrated 
quantum magnets, see, e.g., Refs. [H, 0j| |2o[ |H, [Tslloii 
In particular, in Ref. [l8j it has been demonstrated that 
the CCM GS energy for the s = 1/2 isotropic Heisenberg 
KAFM is close to best available DMRG results : 14 ' 19 

The coupled cluster method (CCM). For the sake of 
brevity we illustrate here only some relevant features of 
the CCM. At that we follow strictly the lines given in 
Ref. Qj|, where the CCM was applied to the isotropic 
spin-s Heisenberg KAFM. For more general information 
on the methodology of the CCM, see, e.g., Refs. l6ll - [67l 
We first mention that the CCM approach yields results 
directly in the thermodynamic limit N —> oo, where N 
is the number of lattice sites. 

First we choose a normalized reference state |<f>) that 
is typically a classical GS of the model. From a quasi- 
classical point of view the coplanar x -\/3 and <7 = 0 
states are favored candidates among the massively de¬ 
generate manifold of classical ground states (see, e.g., 
Refs. 0, il). Consequently, we use both states as ref¬ 
erence states. Then we perform a rotation of the local 
axes of each of the spins such that all spins in the ref¬ 
erence state align along the negative z axis!*! In this 
new set of local spin coordinates we define a complete 
set of mutually commuting multispin creation operators 
Cf = (Gy )t related to this reference state: 

1$) = U-4 •••); C+ = s+ , s+s+ , S+S+S+ , .. . , (2) 

where = s n + I 11 Eq. (j2j) the components of the 

spin operators are defined in the local rotated coordi¬ 
nate frames, and the indices n, m, fc,... denote arbitrary 
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lattice sites, where each site index in each configuration 
index / in Eq. ([2]) can be repeated up to a maximum of 
2s times. With the set {|$), Cf } thus defined, the CCM 
parametrizations of the ket and bra GS eigenvectors |\t f ) 
and ('I'| of the spin system are given by 

|vh) = e 5 |d>) , S = £a/<7+ (3) 

i* o 

(*\ = (<f>\Se~ s , S = l + £>/C 7 . ( 4 ) 

i* o 

The correlation coefficients, aj and a/, contained in the 
CCM correlation operators, S and S, are determined by 
the CCM ket-state and bra-state equations 

($|C7 e~ s He s \<f>) = 0 ; V/ y 2 0 (5) 

($\Se~ s [H, C'+]e s |$) =0 ; VI ^ 0. (6) 

Equations © and © are fully equivalent to the GS 
Schrodinger equations for the ket and bra states. Each 
ket-state or bra-state equation belongs to a certain con¬ 
figuration index /, i.e., it corresponds to a certain set 
(configuration) of lattice sites n, m, k,... , as in Eq. ©. 
Using the Schrodinger equation, H |T) = Eo|T), we can 
now write the GS energy as Eg = (Q\e~ s He s \Q). The 
magnetic order parameter (sublattice magnetization) is 
given by M = —JT = i('I , |sf W; where s\ is expressed 
in the transformed coordinate system, and N{—> oo) is 
the number of lattice sites. 

For the quantum many-body system under considera¬ 
tion we have to use an appropriate approximation scheme 
in order to truncate the expansions of S and S in Eqs. © 
and ©• For that we use the well established SUBn-n 
approximation scheme, cf., e.g., Refs. [HI, [H, HO, HO, lifj - 
16(1 In the SUBro-n scheme we include no more than n 
spin flips spanning a range of no more than n contigu¬ 
ous lattice sites*^ Using an efficient parallelized CCM 
code^ we are able to solve the CCM equations up to 
SUB10-10 for s = 1/2 and up to SUB8-8 for s > 1/2. 
We have calculated the GS energy per spin eg = E 0 /N 
and the magnetic order parameter M for spin quantum 
numbers s = 1/2,1,..., 9/2, 5. The maximum number 
of ket-state equations which we have to take into ac¬ 
count is 416193 for s = 5. Following Ref. [TO we ex¬ 
trapolate the ‘raw’ SUBn-n data to the limit n —> oo 
using n = 4, 5,..., 10 (n = 4, 5,..., 8) for s = 1/2 
(s > 1/2). For that we use the well-tested extrapolation 
ansatz o 12 ' 18 ! 20 ' 48- — eo(n) = ao + ai(l/n) 2 + a 2 (l/n) 4 and 
M(n) = b a + 6i(l/n) 4 / 2 + 6 2 (l/n) 3 / 2 . 

Results. As already mentioned in the introduction, we 
want to calculate the GS s — A phase diagram. Such 
a phase diagram has been very recently presented by 
Chernyshev and Zhitomirsky^ using nonlinear spin-wave 
and real-space perturbation theories. To have the neces¬ 
sary information available for the comparison of Cherny¬ 
shev’s and Zhitomirsky’s results with our CCM results re¬ 
ported below, let us first briefly report the main findings 
of Ref. HJ. It is well know n 2 i 4 i 6 i 18 that in the large-s limit 



FIG. 1. The extrapolated CCM GS energy per spin eo| 7 i->oo 
using the \/3x %/3 reference state (left) and the q = 0 reference 
state (right) as a function of the anisotropy parameter A for 
spin quantum numbers s = 1/2, 1, 3/2,..., 9/2, 5. 
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FIG. 2. The extrapolated CCM GS sublattice magnetization 
Mln-joo using the \/3x y/3 reference state (left) and the q = 0 
reference state (right) as a function of the anisotropy param¬ 
eter A for spin quantum numbers s = 1, 3/2,..., 9/2, 5. Note 
that for s = 1/2 no data are shown, since there is no magnetic 
LRO in the entire region 0 < A < 1. For s = 1 there are finite 
regions of disorder A* < A < 1 for both reference states. 


for A = 1 quantum fluctuations select the x v/3 state. 
In Ref. H3 it was found that this GS selection is preserved 
for weak easy-plane anisotropy, i.e. for A c < A < 1. 
However, the main and unexpected result of Ref. [24] is the 
selection of the q = 0 GS for smaller values of A down to 
the XY point, i.e. for 0 < A < A c . This finding is con¬ 
trary to the selection trend by thermal fluctuations for 
the classical KAFM, where for A = 1 and for A = 0 the 
v/3 x v/3 state is asymptotically selectedi 2 ' 5 ' 71- — Hence, 
for the XY antiferromagnet on the kagome lattice we 
are faced with an example that quantum and thermal 
fluctuations may act very differently. The term in the 
nonlinear spin-wave theory (SWT) responsible for the 
GS selection is of order 0(1), and, therefore, the critical 
anisotropy A c = 0.72235 is found to be independent of 
the spin quantum number s within this approximation* 24 
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FIG. 3. Difference <5e = 0 — e /sx\/3 of the extrapolated 

GS energies of the \/3 x \/3 and the q = 0 states as a function 
of the anisotropy parameter A for spin quantum numbers 
s = 1, 3/2, 5/2, 7/2 and 5. For comparison we also show the 
corresponding large-s results of Ref.)2J (labeled by ’NLSWT’) 
obtained by non-linear SWT. The inset shows an enlarged 
scale of that region of A, where Se changes its sign. 


The real-space perturbation theory provides insight in 
the mechanism of the quantum selection of the ground 
state: Some relevant seventh-order processes change their 
sign as varying A^ 4 The magnetic order parameter cal¬ 
culated in harmonic approximation shows a clear trend 
to magnetic LRO as lowering A. However, this trend 
is substantially overestimated by SWT, since already for 
the extreme quantum case, s = 1/2, the GS is found to 
be magnetically long-range ordered for A < 0.95. On the 
other hand, for the isotropic Heisenberg point it is well- 
known that the linear SWT yields a vanishing order pa¬ 
rameter for any value of si^i Both findings for the order 
parameter are in contradiction to recent results obtained 
by large-scale numerical approaches . 18 ’ 25 

We now discuss our CCM results. In Fig. [Hand Fig. [2] 
we show the GS energy per spin eg and the magnetic 
order parameter M extrapolated to n —> oo for both the 
x/3 x v3 and the q = 0 reference states as a function of 
A for spin quantum numbers up to s = 5. In general, the 
x/3 x -\/3 and the q = 0 cases behave very similar. The 
GS energy for s = 1/2 is almost linearly growing with 
decreasing of A. For larger s the eo(A)-curve noticeably 
deviates from linearity, particularly near A = 1. This 
trend that the influence of the anisotropy A becomes 
exceedingly large near the isotropic Heisenberg limit at 
A = 1 is more pronounced for the order parameter, see 
Fig. [2] For spin quantum numbers s > 3/2 there is a 
drastic downturn in the M (A)-curve as approaching A = 
1 and there is a strong increase in the slope (dM/dA)| a=i 
with growing s. Note that a special behavior for A —> 1 
is also present within the spin wave approach, where the 
1/s corrections diverge for A —> 1 — The cases s = 1/2 



FIG. 4. GS s — A phase diagram of the XXZ spin-s kagome 
antiferromagnet. The open circles connected by the black 
solid line show the critical anisotropy A c (s), where the GS 
selection changes between the %/3 x \/3 state and q = 0 state. 
The upper red line connecting the open squares is the fic¬ 
tional transition line between GS disorder and GS LRO, which 
is obtained by considering the spin quantum number s as a 
continuous quantity. The blue point indicates the transition 
point A* = 0.818 between the disordered ground state and 
the ground state with \/3 x %/3 LRO for s = 1 obtained di¬ 
rectly from the M( A) curve shown in Fig. [2] 


and s = 1 are different from the cases s > 1. For s = 
1/2 our CCM approach leads to a disordered GS in the 
entire region of the anisotropy parameter 0 < A < 1 in 
accordance with recent DMRG calculations.— For s = 1 
we find a finite region of disorder, A* < A < 1, for both 
reference states, where A* = 0.818 (A* = 0.945) for the 
x/3 x >/3 (q = 0) reference state. Note, however, that 
for anisotropies around A* the CCM GS energy for the 
a/ 3 x a/ 3 reference state is lower than that for the q = 0 
reference state, see below. We mention also that a table 
of values of the GS energy and the order parameter for 
A = 1 and for spin quantum numbers up to s = 3 can 
be found in Ref. Il8l . We present corresponding values for 
the XY limit (A = 0) in Table Q] of the present paper. 

Next we discuss the GS selection. For that we consider 
the energy difference Se = eg~° — see Fig. [3] The 

main common feature of the curves shown in Fig.[3]is the 
change of the sign of Se, i.e., in accord with Ref. |24| we 
find a change in the GS selection from the -\/3 x /3 state 
to the g = 0 state at a critical value A c as varying the 
anisotropy from A = 1 to A = 0. The magnitude of 
Se is small, in particular for smaller A. For larger s it 
agrees very well with the SWT data of Ref. [24] up to 
A ~ 0.9. The stronger deviation for A close to one can 
be attributed to the divergence of the 1/s corrections for 
A ->■ 1.— 

As already mentioned above, the large-s spin-wave 
approach^ yields a critical value A c that is independent 
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TABLE I. Extrapolated CCM results for GS energy per spin, 
eo| n ->oo, an d the GS sublattice magnetization, G , using 

the \/3 x \/3 and the q — 0 reference states for A = 0 (XY 
model). Note that corresponding tables for A = 1 (isotropic 
Heisenberg model) can be found in Ref. [3. 




\/3x 

V3 

? = 

0 



eo/s 2 

M/s 

eo/s 2 

M/s 

s = 

1/2 

-1.1896 

< 0 

-1.1968 

< 0 

s = 

1 

-1.0578 

0.8602 

-1.0583 

0.8589 

s = 

3/2 

-1.0347 

0.9402 

-1.0349 

0.9368 

s = 

2 

-1.0252 

0.9570 

-1.0253 

0.9556 

s = 

5/2 

-1.0198 

0.9664 

-1.0199 

0.9656 

s = 

3 

-1.0163 

0.9723 

-1.0164 

0.9719 

s —» 

• OO 

-1 

1 

-1 
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of the spin quantum number s. However, A c certainly 
depends on s. Chernyshev and Zhitomirsky^ suggested 
that A c may increase for smaller spins from the large-s 
value A c = 0.72235 of s — > oo. Our CCM approach yields 
directly A c as a function of s, see the black open circles 
in Fig. 2) Contrary to the conjecture of Chernyshev and 
Zhitomirsky, see Fig. 4b in Ref. [2dl we find that A c in¬ 
creases for larger s. The smallest value, A c = 0.281, is 
found for s = 1. Applying g(x) = a + bx + cx 2 , x = 1/s, 
to extrapolate the CCM data of the critical anisotropy to 
s —> oo yields lim s _ ) . 00 A c = 0.727, where we have used 
s = 3, 7/2,4,9/2, 5 for the extrapolation. This CCM es¬ 
timate of the large-s limit of A c is in excellent agreement 
with the large-s spin-wave result 2i 

To get the full relationship to the GS s — A phase dia¬ 
gram given in Fig. 4 of Ref. [2J we may consider the spin 
quantum number s as a continuous variable. We deter¬ 
mine that (fictional) value of s, for which the GS becomes 
magnetically disordered. For that we fit the data for the 
order parameter M/s as function of s and A by the fit¬ 
ting function /(s, A) = &o(A) — &i(A)s -1 / 2 — & 2 (A)s -1 — 
& 3 (A)s -3 / 2 — & 4 (A)s - 2 , where s = 1/2 is excluded from 
the fit. From /(s c , A) = 0 we obtain the corresponding 
phase boundary s c (A) between magnetically disordered 
and ordered GS phases. This phase boundary is shown 
by the red solid line in Fig. [4j The resulting value of s c 
for A = 1 (A = 0) is s c ~ 1.34 (s c ~ 0.537). 

As already mentioned above for s — 1/2 the GS is 
always magnetically disordered. For s = 1 there are three 
GS phases: the disordered state for 0.818 < A < 1, the 
ordered -\/3 x >/3 state for 0.281 < A < 0.818, and the 
ordered q = 0 state for 0 < A < 0.281. For spin quantum 
numbers s > 1 there are two magnetically ordered GS 
phases, where the phase boundary is given by the black 
solid line in Fig. 2) 


Summary. We summarize our findings by compar¬ 
ing our GS s — A phase diagram (Fig. 21) with that 
of Ref. 24 Most importantly, we get the same trend 
as Ref. [2J, namely the GS selection changes from the 
x/3 x x/3 state to the q = 0 state at a critical A c as vary¬ 
ing the anisotropy parameter A from A = 1 to A = 0. 
The energy difference <5e between the \/3 x \/3 and the 
<7 = 0 state calculated by the CCM, e.g. for s = 5, is in 
excellent agreement with the large-s SWT data of Ref. l24l. 
Moreover, we find also in accordance with Ref. [24] that 
the region of disorder in the GS s — A phase diagram 
grows with growing A. 

Since our CCM approach is not limited to large s, we 
can overcome some limitations of the SWT of Ref. [24|: 
The CCM result for the critical A c depends on s by con¬ 
trast to the s- independent SWT value. For s = 1/2 the 
SWT yields GS magnetic LRO for A < 0.95, whereas the 
CCM (in agreement with DMRG results^) yields always 
a disordered GS. 

Knowing these limitations of the SWT, Chernyshev 
and Zhitomirsky^ proposed a tentative GS s — A phase 
diagram, see Fig. 4b of Ref. HJ. However, our phase dia¬ 
gram determined from the CCM results differs from that 
conjectured by Chernyshev and Zhitomirsky. We find 
that A c increases with increasing s, whereas Cherny¬ 
shev and Zhitomirsky speculate that A c may decrease 
for larger s. Hence our results indicate that the region of 
x/3 x a/ 3 GS LRO is much larger than that proposed in 
Ref. [24]. As a consequence, the (fictional) transition from 
the magnetically disordered GS to a state with magnetic 
LRO obtained by a continuous increase of spin quan¬ 
tum number s starting from the extreme quantum case 
s = 1/2 is always a transition to the state with /3 x y/3 
LRO, while in Ref. [U it is suggested that the transition 
is to the state with q = 0 LRO. Moreover, for s = 3/2 we 
find 1/3 x a/ 3 LRO for 0.525 < A < 1, while in Ref. [24| 
q = 0 LRO is suggested for the entire region 0 < A < 1. 

We conclude, that the interplay of frustration, quan¬ 
tum fluctuations and anisotropy leads to a rich ground- 
state phase diagram of the XXZ spin-s KAFM. Bear¬ 
ing in mind the numerous investigations of the isotropic 
Heisenberg KAFM, see, e.g., Refs. Iffl and references 
therein, the anisotropic model provides a challenging 
playground to apply the toolbox of frustrated quantum 
magnetism on this so far little investigated problem. 
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